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Abstract — Humidity  over  the  ocean  tends  to  be  higher  than 
over  land,  and  long-range  links  at  frequencies  above  10  GHz  may 
experience  large  propagation  loss  due  to  atmospheric  absorption. 
When  high  link  reliability  is  required,  the  log-normal  assumption 
used  to  estimate  power  margins  from  the  loss  mean  and  standard 
deviation  is  not  accurate.  To  solve  these  problems,  a  set  of 
atmospheric  models  has  been  developed  from  radiosonde  profiles 
collected  along  the  Atlantic  coast  of  the  US.  The  large  quantity  of 
radiosonde  data  available  was  used  to  identify  radiosonde  profiles 
that  occur  with  specified  probabilities  at  different  times  of  year 
and  locations.  The  result  was  used  to  estimate  high-reliability 
SHF/EHF  air-to-surface  radio  link  performance  in  a  maritime 
environment. 

I.  Introduction 

Air-to-surface  radio  links  differ  from  typical  satellite  com¬ 
munications  links  in  that  the  path  elevation  angles  are  lower, 
often  2°  —  5°,  and  have  a  longer  extent  through  the  atmosphere. 
When  signal  frequencies  exceed  10  GHz,  propagation  loss 
due  to  absorption  by  atmospheric  gases  becomes  a  significant 
link  budget  term.  Losses  over  the  ocean  are  even  higher 
than  over  land,  due  to  increased  humidity.  Finally,  when 
high  reliability  is  desired,  the  log-normal  assumption  used 
to  estimate  required  power  margins  from  the  loss  mean  and 
standard  deviation  is  not  accurate.  A  set  of  atmospheric  loss 
models  has  been  developed  from  radiosonde  profiles  collected 
along  the  Atlantic  coast  of  the  United  States,  to  accurately 
estimate  high-reliability  SHF/EHF  air-to-surface  radio  link 
performance  in  a  maritime  environment.  Data  was  collected 
at  three  locations  approximately  spanning  the  coast:  Chatham 
MA,  Wallops  Island  VA,  and  Key  West,  FL.  Weather  station 
locations  and  identifiers  are  listed  in  Table  I.  The  large  quantity 
of  historical  radiosonde  data  available  for  those  locations  was 
used  to  identify  radiosonde  profiles  that  occur  with  specified 
probabilities.  Those  profiles  were  used  to  predict  required  link 
performance  to  achieve  high  reliability  at  different  locations 
and  times  of  year. 

II.  Data  Acquisition 

Radiosonde  balloons  are  launched  daily  at  selected  loca¬ 
tions,  and  measure  temperature,  dew  point  temperature,  and  air 
pressure  as  they  ascend.  Radiosondes  are  ordinarily  launched 
twice  a  day,  nominally  at  midnight  and  noon  UTC,  (7  AM  and 
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TABLE  I.  Weather  Stations 


Location 

WMO 

Number 

Station 

ID 

Latd. 

Long. 

Alt. 

Chatham, 

Massachusetts 

74494 

CHH 

41.67° 

69.97° 

16  m 

Wallops  Island, 
Virginia 

72402 

WAL 

37.85° 

75.48° 

3  m 

Key  West, 
Florida 

72201 

EYW 

24.55° 

81.75° 

1  m 

7  PM  EST,  8  AM  and  8  PM  EDT).  Actual  launch  time  is  usu¬ 
ally  less  than  an  hour  before  nominal  launch  time.  Radiosonde 
data  is  available  on  the  internet  at  http://esrl.noaa.gov/roabs, 
for  dates  listed  in  [1].  The  data  format  is  described  in  [2],  [3], 
[4].  The  models  presented  here  are  based  on  available  files 
from  2004  -  2013.  Text  files  containing  radiosonde  data  for 
one  month  of  each  year  were  downloaded  and  converted  to 
Matlab  matrices,  then  processed  as  follows: 

1)  Unusable  or  undesirable  lines  were  deleted,  e.g.  lines 
with  missing  measurements,  or  multiple  copies  of  a 
single  line. 

2)  Altitude  measurements  were  converted  from  geopo¬ 
tential1  altitude  to  geometric  altitude  [5]. 

3)  Water  vapor  partial  pressure  was  calculated  at  each 
altitude  from  temperature,  dew  point  temperature,  and 
air  pressure  [6]. 

4)  Water  vapor  density  was  calculated  from  water  vapor 
partial  pressure  and  temperature  [7]. 

5)  Data  was  interpolated  to  a  set  of  consistent  and  more 
closely  spaced  altitudes 

III.  Propagation  Loss  Calculations 

The  analysis  used  assumes  weather  data  are  constant  inside 
nested  earth-centered  spherical  shells,  with  shell  boundaries 
at  altitudes  recommended  in  [8].  In  step  (5)  above,  data  was 
interpolated  to  the  mid-point  of  each  shell.  Using  weather  data, 
the  index  of  refraction,  or  equivalently  complex  refractivity 
can  be  found  in  each  atmosphere  shell.  The  real  part  of  the 
index  of  refraction  determines  how  the  radio  path  bends  at 
shell  boundaries,  and  is  easily  found  to  the  necessary  accuracy 
[9].  The  imaginary  part  of  refractivity  determines  atmospheric 
absorption,  and  was  more  carefully  approximated  using  Liebe’s 
approximation  [10],  which  is  good  for  frequencies  below  100 
GHz.  It  accounts  for  absorption  by  water  vapor  and  oxygen. 

1  Geopotential  altitude  is  defined  by  a  surface  of  constant  gravitational 
potential  energy,  and  does  not  equal  geometric  altitude.  At  20  km  altitude, 
the  difference  is  120  m  [5]. 


Water  vapor  has  a  resonant  line  at  22.2  GHz  and  dominates  loss 
in  the  lower  end  of  the  band,  oxygen  has  a  series  of  resonant 
lines  around  60  GHz,  and  makes  a  significant  loss  contribution 
at  higher  frequencies.  The  version  of  Liebe’s  approximation 
used  is  in  [8].  Path  length  was  found  using  a  ray-tracing  model 
[8].  Starting  at  ground  level,  the  angle  at  which  an  uplink  path 
enters  a  new  shell  was  found  from  the  angle  at  which  it  left 
the  lower  shell  and  the  shell  indices  of  refraction.  The  path 
was  traced  upward  until  it  reached  the  desired  altitude,  and 
the  ground  range  was  found  from  the  central  angle  subtended 
by  the  path.  Gas  absorption  along  the  path  was  found  as 
the  summed  product  of  path  length  through  each  shell,  and 
the  specific  attenuation  (dB/km)  in  each  shell  [8].  Specific 
attenuation  is  a  simple  function  of  the  imaginary  refractivity 
[8],  [10]. 


a.  Chatham,  MA 


IV.  Atmospheric  Characterization 


Ten  years  of  data  were  collected  and  analyzed  for  each 
location.  Figure  1  shows  plots  of  parameter  means  and  stan¬ 
dard  deviations  versus  altitude  at  Chatham  and  Key  West 
for  July.  The  region  between  p  ±  a  is  colored  light  green. 
ITU  mid-latitude  standard  atmospheres  [11]  were  also  plotted. 
The  Chatham  July  mean  temperature  falls  close  to  the  ITU 
summer  model,  slightly  below  at  the  surface,  and  slightly 
above  elsewhere.  The  Key  West  July  mean  temperature  is 
more  than  one  standard  deviation  above  the  ITU  summer 
model  at  all  altitudes.  Air  pressure  variance  at  both  locations 
is  small  compared  to  the  variation  with  altitude,  and  the  green 
region  is  mostly  concealed  behind  the  ITU  summer  model. 
Mean  water  vapor  partial  pressure  exceeds  the  ITU  Summer 
model  by  less  than  one  standard  deviation  at  Chatham,  and 
by  several  times  the  standard  deviation  at  Key  West.  Wallops 
Island  results  fall  between  Chatham  and  Key  West.  Converting 
the  large  quantity  of  radiosonde  traces  to  propagation  loss, 
then  collecting  loss  statistics  is  a  time  consuming  process. 
Furthermore,  atmosphere  loss  margins  for  high  reliability  links 
calculated  from  those  statistics  under  the  assumption  of  log- 
normally  distributed  loss  are  inaccurate.  As  an  alternative,  all 
radiosonde  data  collects  for  a  given  month  over  ten  years 
have  been  ranked  according  to  their  approximate  severity  of 
propagation  loss.  Data  files  with  the  selected  rank  were  used 
to  find  loss  margins,  i.e.,  atmosphere  loss  calculated  from  the 
99th  percentile  data  file  was  taken  as  the  margin  required  for 
99%  reliability.  Radiosonde  profiles  were  ranked  according  to 
their  columnar  water  vapor,  defined  as  the  total  mass  of  water 
vapor  in  a  cylinder  with  a  1  m2  face,  extending  vertically  from 
ground  level  to  the  aircraft  altitude.  The  columnar  water  vapor 
is 

_ _ _  Ah  n- 

(1) 


5hipi 

c  =  2^Wokg 


where  Shi  is  the  thickness  of  shell  i,  and  pi  is  water  vapor 
density  ( g/m 3)  inside  shell  i. 

Figure  2  shows  scatter  plots  of  atmoshpere  loss  versus 
columnar  water  vapor  at  two  frequencies,  for  a  path  with 
aircraft  altitude  35,000  feet  (10.7  km),  and  ground  range 
120  nmi  (222  km).  At  this  altitude,  120  nmi  ground  range 
corresponds  to  a  2°  elevation  angle.  Data  is  shown  for  all 
three  locations  and  every  season. 


b.  Key  West,  FL 

Fig.  1.  July  weather  statistics  from  radiosonde  data,  2004-2013. 


nearly  perfect  straight  line.  At  44  GHz,  oxygen  absorption 
is  also  a  factor.  Oxygen  absorption  varies  inversely  with 
temperature,  which  has  little  effect  on  columnar  water  vapor. 
However,  oxygen  loss  is  less  variable  than  water  vapor  loss, 
and  columnar  water  vapor  was  considered  adequate  for  ranking 
weather  profiles  at  44  GHz. 

V.  The  Chatham  Model 

Figure  3  shows  smoothed  histograms  of  columnar  water 
vapor  in  Chatham  for  each  season.  In  Chatham  cool,  dry 
months  display  less  variation  than  warm,  humid  months,  as 
probability  accumulates  near  the  minimum  value  of  zero.  At 
Key  West  the  opposite  effect  is  observed,  and  hot,  wet  months 
show  lower  variability  as  probability  accumulates  near  the 
maximum  value  obtained  when  humidity  is  100%. 

Figure  4  plots  water  vapor  density  versus  altitude  for  some 
profiles  from  Chatham  in  July,  ranked  by  their  columnar  water 
vapor.  Figure  5  is  a  plot  of  atmospheric  loss  versus  ground 
range  at  21.7  GHz,  for  the  radiosonde  profiles  of  Figure  4. 
It  also  shows  loss  for  the  maximum  and  minimum  columnar 
water  vapor  during  the  ten-year  interval.  Similar  results  were 
found  at  44  GHz.  The  44  GHz  median  loss  at  120  nmi  (222 
km)  was  slightly  higher,  but  the  99th  percentile  loss  at  that 
range  was  about  3  dB  lower.  The  spread  was  reduced  at  44 
GHz  because  oxygen  absorption  is  less  variable. 


Water  vapor  absorption  is  the  dominant  loss  mechanism 
near  22  GHz,  and  the  21  GHz  data  in  Figure  2a  forms  a 


Figure  6  shows  annual  variation  of  columnar  water  vapor 
at  Chatham,  for  representative  percentiles.  Figure  7  shows  gas 


a.  20.7  GHz 


Columnar  Water  Vapor  (kg) 

b.  44.5  GHz 

Fig.  2.  Gas  loss  versus  columnar  water  vapor.  Ground  range  120  nmi  (222 
km),  aircraft  altitude  35,000  feet  (10.7  km). 


loss  to  an  aircraft  at  35,000  feet  (10.7  km)  altitude  with  120 
nmi  (222  km)  ground  range,  using  the  radiosonde  profiles 
corresponding  to  the  percentiles  in  Figure  6.  The  21  GHz 
median  loss  varies  by  about  12  dB  over  a  year,  and  99th 
percentile  loss  varies  by  about  14  dB.  Loss  curves  for  the 
44  GHz  signal  sometimes  cross  each  other  because  files  were 
ranked  by  columnar  water  vapor,  which  does  not  account  for 
oxygen  absorption.  Oxygen  absorption  varies  inversely  with 
temperature,  and  may  be  high  on  cool  dry  days  when  columnar 
water  vapor  is  low.  The  44  GHz  annual  variation  is  8-9  dB  for 
the  median,  and  about  12  dB  for  the  99th  percentile. 

VI.  Other  Locations 

Weather  data  was  analyzed  similarly  for  Wallops  Island  and 
Key  West.  Figure  8  shows  annual  variation  in  columnar  water 
vapor  at  Wallops  Island.  Figure  9  shows  gas  loss  at  120  nmi 
(222  km)  ground  range  to  an  aircraft  at  35,000  feet  (10.7  km) 
altitude,  using  the  same  radiosonde  profiles  as  Figure  8.  Figure 
10  shows  columnar  water  vapor  percentiles  by  month  at  Key 


Fig.  3.  Columnar  water  vapor  histograms  from  Chatham  radiosonde  data, 
2004-2013. 


West,  and  Figure  11  shows  corresponding  gas  loss  percentiles. 
The  annual  variation  in  median  gas  loss  is  14  dB  at  21  GHz, 
and  9  dB  at  44  GHz.  The  annual  variation  to  the  99th  percentile 
line  is  less,  10  dB  and  7  dB  respectively. 

Figure  12  compares  propagation  losses  at  all  locations  in 
July,  the  worst  or  near-worst  month  at  each  location,  versus 
reliability.  Reliability  is  the  probability  that  observed  loss  will 
be  no  greater  than  the  value  indicated.  For  the  data  presented 
here,  loss  for  a  given  reliability  was  estimated  to  equal  the 
measured  loss  percentile.  Curves  are  more  variable  at  higher 
reliability  because  those  estimates  were  based  on  events  in  the 
tails  of  the  distributions,  where  there  was  less  data. 

Comparing  both  plots  in  Figure  12,  loss  increases  as  the 
location  moves  south  to  hotter,  more  humid  climates  but 
atmospheric  loss  variation  decreases.  The  difference  between 
locations  is  most  pronounced  at  low  reliability,  becoming 
smaller  as  reliability  increases.  Wallops  Island  and  Key  West 


Water  Vapor  Density  (g/m3) 


Fig.  4.  Water  vapor  density  profiles  from  radiosonde  data  with  indicated 
columnar  water  vapor  rank,  measured  over  altitudes  0  -  35,000  feet  (10.7 
km),  Chatham  2004-2013. 


Ground  Range  (nmi) 


Fig.  5.  Atmospheric  loss  versus  ground  range  for  indicated  columnar  water 
vapor  rank.  Frequency  21.2  GHz,  aircraft  altitude  35,000  feet  (10.7  km), 
Chatham,  July,  2004-2013. 

have  almost  identical  downlink  losses  at  99%  reliability, and 
Chatham  loss  is  only  2-3  dB  lower.  In  comparison  to  the  results 
shown  here,  the  worst  case  loss  at  X-band  (8.5  GHz)  is  1.5  dB. 
This  figure  demonstrates  once  again  that  there  is  a  maximum 
vapor  pressure  which  as  approached  in  hot,  wet  climates, 
limiting  the  maximum  gas  absorption  that  may  be  observed. 
An  increase  in  humidity  beyond  the  maximum  results  in  clouds 
or  fog,  loss  mechanisms  that  are  not  treated  in  this  paper. 

VII.  Conclusion 

The  availability  of  weather  service  data  on  the  internet  is 
a  valuable  resource  for  predicting  performance  of  radio  links 
above  10  GHz.  Radiosonde  profiles  can  be  used  to  predict 
atmospheric  propagation  for  a  particular  time  and  place.  His¬ 
torical  collections  of  data  allow  us  to  predict  power  margins  re¬ 
quired  to  compensate  for  infrequently  occurring  conditions,  by 
ranking  past  radiosonde  profiles  and  identifying  those  causing 
high  loss.  As  the  quantity  of  historical  data  grows,  increasingly 
improbable  conditions  can  be  modeled.  This  information  can 
also  be  combined  with  weather  prediction  to  plan  for  upcoming 
severe  propagation  weather  events  for  operational  systems.  As 
an  example,  available  data  from  three  weather  stations  along 
the  east  coast  of  the  US  were  used  to  develop  models  that 
predict  SHF/EHF  air-to-ground  link  performance  in  a  maritime 
environment.  It  was  found  that  although  there  is  significant 
variation  in  median  performance  from  north  to  south,  the 
margins  required  to  achieve  greater  than  90%  reliability  vary 
by  only  few  dB  between  Massachusetts  and  the  southern  tip 
of  Florida. 


Fig.  6.  Columnar  water  vapor  percentiles  by  month,  Chatham,  2004-2013. 
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Fig.  7.  Atmospheric  loss  for  indicated  columnar  water  vapor  percentiles 
by  month,  altitude  35,000  feet  (10.7  km),  ground  range  120  nmi  (222  km), 
Chatham  radiosonde  data  2004-2013. 


Fig.  8.  Columnar  vapor  percentiles  by  month,  Wallops  Island,  2004-2013. 


a.  20.7  GHz 


Fig.  9.  Atmospheric  loss  for  indicated  columnar  water  vapor  percentiles 
by  month,  altitude  35,000  feet  (10.7  km),  ground  range  120  nmi  (222  km), 
Wallops  Island  radiosonde  data  2004-2013. 


Fig.  10.  Columnar  vapor  percentiles  by  month,  Key  West,  2004-2013. 
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Fig.  11.  Atmospheric  loss  for  indicated  columnar  water  vapor  percentiles  by 
month,  altitude  35,000  feet  (10.7  km),  ground  range  120  nmi  (222  km),  Key 
West  radiosonde  data  2004-2013. 
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Fig.  12.  July  gas  loss  margin  versus  reliability,  parameterized  by  ground 
range.  Aircraft  altitude  35,000  feet  (10.7  km). 


